- Author:
- nima <nafs080@aucklanduni.ac.nz>
- Date:
- 2020-07-17 14:20:26+12:00
- Desc:
- test-SEDML
- Permanent Source URI:
- https://models.fieldml.org/workspace/572/rawfile/75b8e504afd4e0a405f0012421e3dfce362a77aa/New - Composite Model(27.02.2019).cellml
<?xml version='1.0' encoding='UTF-8'?>
<model name="Nima_2018" xmlns="http://www.cellml.org/cellml/1.0#" xmlns:cellml="http://www.cellml.org/cellml/1.0#">
<!-- Unit and value for K_s (GLUT2) is fixed
NHE3 and NBC were fixed also sign of J_NHE3 changed from + to -
decrease D_P_Na from 1e-8 t0 1e-9
JENaC was fixed (-to+) and the sign for J_Enac and I_Enac changed from - to +
CFTR was fixed (-to+) nad the sign for J_Enac and I_Enac changed from - to + -->
<units name="per_second">
<unit exponent="-1" units="second"/>
</units>
<units name="per_second2">
<unit exponent="-1" units="second"/>
</units>
<!-- def unit per_second5 as
unit second {expo: -5};
enddef;-->
<units name="volt_per_second">
<unit units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_meter">
<unit exponent="-1" units="meter"/>
</units>
<units name="per_volt_per_second">
<unit exponent="-1" units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_volt4_per_second">
<unit exponent="-4" units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_volt3_per_second">
<unit exponent="-3" units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_volt2_per_second">
<unit exponent="-2" units="volt"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_mol">
<unit exponent="-1" units="mole"/>
</units>
<units name="umol">
<unit prefix="micro" units="mole"/>
</units>
<units name="per_umol">
<unit exponent="-1" units="umol"/>
</units>
<units name="M">
<unit units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="per_M">
<unit exponent="-1" units="M"/>
</units>
<units name="m2">
<unit exponent="2" units="metre"/>
</units>
<units name="cm">
<unit prefix="centi" units="metre"/>
</units>
<units name="um2">
<unit exponent="2" prefix="micro" units="metre"/>
</units>
<units name="um">
<unit prefix="micro" units="meter"/>
</units>
<units name="um3">
<unit exponent="3" prefix="micro" units="metre"/>
</units>
<units name="m3">
<unit exponent="3" units="metre"/>
</units>
<!-- def unit litre as
unit m3 {pref: milli};
enddef;-->
<units name="mol_per_um2">
<unit units="mole"/>
<unit exponent="-1" units="um2"/>
</units>
<units name="mM">
<unit prefix="milli" units="mole"/>
<unit exponent="-1" units="litre"/>
</units>
<units name="mV">
<unit prefix="milli" units="volt"/>
</units>
<units name="J_per_K_per_mol">
<unit exponent="-1" units="mole"/>
<unit units="joule"/>
<unit exponent="-1" units="kelvin"/>
</units>
<units name="pA_per_cm2">
<unit prefix="pico" units="ampere"/>
<unit exponent="-2" prefix="centi" units="metre"/>
</units>
<units name="uA">
<unit prefix="micro" units="ampere"/>
</units>
<units name="uA_per_um2">
<unit units="uA"/>
<unit exponent="-1" units="um2"/>
</units>
<units name="uA_per_A">
<unit units="uA"/>
<unit exponent="-1" units="ampere"/>
</units>
<units name="umol_per_mol">
<unit units="umol"/>
<unit exponent="-1" units="mole"/>
</units>
<units name="C_per_mol">
<unit exponent="-1" units="mole"/>
<unit units="coulomb"/>
</units>
<units name="per_M_per_second">
<unit exponent="-1" units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_M2_per_second">
<unit exponent="-2" units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_M4_per_second2">
<unit exponent="-4" units="M"/>
<unit exponent="-2" units="second"/>
</units>
<units name="per_M4_per_second">
<unit exponent="-4" units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="uM">
<unit prefix="micro" units="M"/>
</units>
<units name="M2">
<unit exponent="2" units="M"/>
</units>
<units name="M3">
<unit exponent="3" units="M"/>
</units>
<units name="M3_per_second">
<unit exponent="3" units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="per_M3_per_second5">
<unit exponent="-3" units="M"/>
<unit exponent="-5" units="second"/>
</units>
<units name="M_per_second">
<unit units="M"/>
<unit exponent="-1" units="second"/>
</units>
<units name="cm2">
<unit exponent="2" prefix="centi" units="metre"/>
</units>
<units name="cm2_per_s">
<unit units="cm2"/>
<unit exponent="-1" units="second"/>
</units>
<units name="mole_per_m2">
<unit units="mole"/>
<unit exponent="-1" units="m2"/>
</units>
<!-- def unit uM_per_cm2_per_s as
unit uM;
unit cm2 {expo: -1};
unit second {expo: -1};
enddef;-->
<units name="Farad">
<unit exponent="-1" units="volt"/>
<unit units="coulomb"/>
</units>
<units name="uF">
<unit prefix="micro" units="Farad"/>
</units>
<units name="uF_per_cm2">
<unit prefix="micro" units="Farad"/>
<unit exponent="-1" units="cm2"/>
</units>
<units name="umol_per_cm2_per_s">
<unit prefix="micro" units="mole"/>
<unit exponent="-1" units="cm2"/>
<unit exponent="-1" units="second"/>
</units>
<units name="umol_per_s">
<unit prefix="micro" units="mole"/>
<unit exponent="-1" units="second"/>
</units>
<units name="uA_per_cm2">
<unit units="uA"/>
<unit exponent="-2" prefix="centi" units="metre"/>
</units>
<units name="m_per_s">
<unit units="meter"/>
<unit exponent="-1" units="second"/>
</units>
<units name="m3_mole_per_litre_per_umol">
<unit units="m3"/>
<unit units="mole"/>
<unit exponent="-1" units="litre"/>
<unit exponent="-1" units="umol"/>
</units>
<units name="litre_umol_per_mol_per_m3">
<unit units="litre"/>
<unit units="umol"/>
<unit exponent="-1" units="mole"/>
<unit exponent="-1" units="m3"/>
</units>
<units name="S">
<unit units="siemens"/>
</units>
<units name="S_per_m2">
<unit units="siemens"/>
<unit exponent="-1" units="m2"/>
</units>
<units name="per_m2">
<unit exponent="-1" units="m2"/>
</units>
<units name="cm_per_s">
<unit units="cm"/>
<unit exponent="-1" units="second"/>
</units>
<component name="parameters">
<variable name="time" public_interface="out" units="second"/>
<variable initial_value="-1" name="z_Cl" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="z_Na" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="z_k" public_interface="out" units="dimensionless"/>
<variable initial_value="96485" name="F" public_interface="out" units="C_per_mol"/>
<variable initial_value="8.314" name="R" public_interface="out" units="J_per_K_per_mol"/>
<variable initial_value="310" name="T" public_interface="out" units="kelvin"/>
<variable initial_value="1e-5" name="capacitance" public_interface="out" units="uF"/>
<variable initial_value="8.4e-10" name="A_Apical" public_interface="out" units="m2"/>
<variable initial_value="5e-10" name="A_Basol" public_interface="out" units="m2"/>
<variable initial_value="2e-15" name="v_cell" public_interface="out" units="m3"/>
<variable initial_value="0" name="P_Na" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_k" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_Cl" public_interface="out" units="m_per_s"/>
<variable initial_value="1e-8" name="P_P_Na" public_interface="out" units="m_per_s"/>
<variable initial_value="5e-8" name="P_P_k" public_interface="out" units="m_per_s"/>
<variable initial_value="2.2e-8" name="P_P_Cl" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_Bl_Na" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_Bl_k" public_interface="out" units="m_per_s"/>
<variable initial_value="0" name="P_Bl_Cl" public_interface="out" units="m_per_s"/>
<variable initial_value="14000" name="k0_12" public_interface="out" units="per_M2_per_second"/>
<variable initial_value="300" name="k0_21" public_interface="out" units="per_second"/>
<variable initial_value="15000" name="k0_23" public_interface="out" units="per_M_per_second"/>
<variable initial_value="20" name="k0_32" public_interface="out" units="per_second"/>
<variable initial_value="0.1" name="k0_25" public_interface="out" units="per_second"/>
<variable initial_value="0.01" name="k0_52" public_interface="out" units="per_second"/>
<variable initial_value="50" name="k0_34" public_interface="out" units="per_second"/>
<variable initial_value="50" name="k0_43" public_interface="out" units="per_second"/>
<variable initial_value="6" name="k0_45" public_interface="out" units="per_second"/>
<variable initial_value="450" name="k0_54" public_interface="out" units="per_M_per_second"/>
<variable initial_value="10" name="k0_56" public_interface="out" units="per_second"/>
<variable initial_value="200" name="k0_65" public_interface="out" units="per_M2_per_second"/>
<variable initial_value="25" name="k0_61" public_interface="out" units="per_second"/>
<variable initial_value="600" name="k0_16" public_interface="out" units="per_second"/>
<variable initial_value="0.7" name="delta" public_interface="out" units="dimensionless"/>
<variable initial_value="0.3" name="alpha_p" public_interface="out" units="dimensionless"/>
<variable initial_value="0" name="alpha_pp" public_interface="out" units="dimensionless"/>
<variable initial_value="2" name="n" public_interface="out" units="dimensionless"/>
<variable initial_value="-2" name="z_c" public_interface="out" units="dimensionless"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="mu" public_interface="out" units="dimensionless"/>
<!--alpha_p = 1{dimensionless}-delta-alpha_pp;-->
<!--C_T = 1e6{umol_per_mol}*N_C/N_Avo;-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>mu</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="Apical_concentrations">
<variable initial_value="0.105" name="Na_m" public_interface="out" units="M"/>
<variable initial_value="0.004" name="K_m" public_interface="out" units="M"/>
<variable initial_value="0.05" name="glucose_m" public_interface="out" units="M"/>
<variable initial_value="0.106" name="Cl_m" public_interface="out" units="M"/>
<!-- var k: per_second {init: 0.2, pub: out};
var glucose_0: M {init: 20e-3, pub: out};
var glucose_f: M {init: 0.1e-3, pub: out};
var time: second {pub: in};-->
<!-- glucose_m = sel
case (time >= 0{second}) and (time < 60{second}):
0.1e-3{M};
case time >= 60{second}:
(glucose_0-glucose_f)*exp(-k*(time-60{second}))+glucose_f;
endsel;-->
</component>
<component name="Apical_voltage">
<variable initial_value="-0.046" name="v_mc" public_interface="out" units="volt"/>
<variable name="time" public_interface="in" units="second"/>
<variable name="capacitance" public_interface="in" units="uF"/>
<variable name="I_D_P_Na" public_interface="in" units="uA"/>
<variable name="I_D_P_K" public_interface="in" units="uA"/>
<variable name="I_D_P_Cl" public_interface="in" units="uA"/>
<variable name="I_NHE3_H" public_interface="in" units="uA"/>
<variable name="I_NHE3_Na" public_interface="in" units="uA"/>
<variable name="I_BK" public_interface="in" units="uA"/>
<!-- var I_AE1_HCO3: uA {pub: in};
var I_AE1_Cl: uA {pub: in};-->
<variable name="I_ENaC" public_interface="in" units="uA"/>
<!--var I_NaCl: uA {pub: in};-->
<variable name="I_Na_SGLT" public_interface="in" units="uA"/>
<variable name="I_CFTR" public_interface="in" units="uA"/>
<!--var I_CFTR_HCO3: uA {pub: in};-->
<!--var I_Gl_SGLT: uA {pub: in};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v_mc</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_Na_SGLT</ci>
</apply>
<ci>I_BK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_CFTR</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_ENaC</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<ci>I_D_P_Na</ci>
<ci>I_D_P_K</ci>
<ci>I_D_P_Cl</ci>
</apply>
</apply>
</apply>
</apply>
<ci>capacitance</ci>
</apply>
</apply>
</math>
</component>
<component name="Cell_concentration">
<variable initial_value="0.047" name="Na_i" public_interface="out" units="M"/>
<variable initial_value="0.003" name="glucose_i" public_interface="out" units="M"/>
<variable initial_value="0.129" name="K_i" public_interface="out" units="M"/>
<variable initial_value="6e-5" name="Ca_i" public_interface="out" units="M"/>
<variable initial_value="0.057" name="Cl_i" public_interface="out" units="M"/>
<variable name="HCO3_i" public_interface="out" units="M"/>
<variable name="H_int" public_interface="out" units="M"/>
<variable initial_value="0.003" name="HCO3_m" public_interface="out" units="M"/>
<variable initial_value="0.041" name="HCO3_s" public_interface="out" units="M"/>
<variable name="H_Ext" public_interface="out" units="M"/>
<variable name="H_s" public_interface="out" units="M"/>
<variable name="v_cell" public_interface="in" units="m3"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="J_Na_SGLT" public_interface="in" units="umol_per_s"/>
<variable name="J_Gl_SGLT" public_interface="in" units="umol_per_s"/>
<variable name="J_Na_NaK" public_interface="in" units="umol_per_s"/>
<variable name="J_K_NaK" public_interface="in" units="umol_per_s"/>
<variable name="J_NBC_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_NBC_HCO3" public_interface="in" units="umol_per_s"/>
<!--var J_NaCl: umol_per_s {pub: in};-->
<variable name="J_NHE3_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_NHE3_H" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Bl_Na" public_interface="in" units="umol_per_s"/>
<variable name="J_GLUT" public_interface="in" units="umol_per_s"/>
<variable name="J_A_GLUT" public_interface="in" units="umol_per_s"/>
<variable name="J_D_K" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Bl_K" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Cl" public_interface="in" units="umol_per_s"/>
<variable name="J_D_Bl_Cl" public_interface="in" units="umol_per_s"/>
<variable name="J_AE1_Cl" public_interface="in" units="umol_per_s"/>
<variable name="J_AE1_HCO3" public_interface="in" units="umol_per_s"/>
<variable name="J_CFTR" public_interface="in" units="umol_per_s"/>
<variable name="J_CLC2" public_interface="in" units="umol_per_s"/>
<variable name="J_CFTR_HCO3" public_interface="in" units="umol_per_s"/>
<variable name="J_ENaC" public_interface="in" units="umol_per_s"/>
<variable name="J_NaKCC" public_interface="in" units="umol_per_s"/>
<variable name="J_IK" public_interface="in" units="umol_per_s"/>
<variable name="J_BK" public_interface="in" units="umol_per_s"/>
<variable name="J_D_P_Na" public_interface="in" units="umol_per_s"/>
<!-- var J_buf_C: M_per_second {pub: in};
var J_buf_A: M_per_second {pub: in};-->
<variable name="time" public_interface="in" units="second"/>
<variable initial_value="7.2" name="pH_int" public_interface="out" units="per_second"/>
<variable initial_value="7.2" name="pH_Ext" public_interface="out" units="per_second"/>
<variable initial_value="7.2" name="pH_s" public_interface="out" units="per_second"/>
<variable initial_value="1680000" name="rho" public_interface="out" units="per_meter"/>
<variable name="beta" public_interface="out" units="M"/>
<variable initial_value="6.3" name="p_K" public_interface="out" units="dimensionless"/>
<variable initial_value="40" name="P_CO2" public_interface="out" units="dimensionless"/>
<variable initial_value="1e-5" name="s" public_interface="out" units="dimensionless"/>
<variable initial_value="1" name="theta_1" units="dimensionless"/>
<variable initial_value="1" name="theta_2" units="dimensionless"/>
<variable initial_value="1" name="theta_3" units="dimensionless"/>
<variable initial_value="1" name="theta_4" units="dimensionless"/>
<variable initial_value="1" name="theta_5" units="dimensionless"/>
<variable initial_value="1" name="theta_6" units="dimensionless"/>
<variable initial_value="1" name="theta_7" units="dimensionless"/>
<variable initial_value="1" name="theta_8" units="dimensionless"/>
<variable initial_value="1" name="theta_9" units="dimensionless"/>
<variable initial_value="1" name="theta_10" units="dimensionless"/>
<variable initial_value="1" name="theta_11" units="dimensionless"/>
<variable initial_value="1" name="theta_12" units="dimensionless"/>
<variable initial_value="1" name="theta_13" units="dimensionless"/>
<variable initial_value="1" name="theta_14" units="dimensionless"/>
<variable initial_value="1" name="theta_15" units="dimensionless"/>
<variable initial_value="1" name="theta_16" units="dimensionless"/>
<variable initial_value="1" name="theta_17" units="dimensionless"/>
<variable initial_value="1" name="theta_18" units="dimensionless"/>
<variable initial_value="1" name="theta_19" units="dimensionless"/>
<variable initial_value="1" name="theta_20" units="dimensionless"/>
<variable initial_value="1" name="theta_21" units="dimensionless"/>
<variable initial_value="1" name="theta_22" units="dimensionless"/>
<variable initial_value="1" name="theta_23" units="dimensionless"/>
<variable initial_value="1" name="theta_24" units="dimensionless"/>
<variable initial_value="1" name="theta_25" units="dimensionless"/>
<variable initial_value="1" name="theta_26" units="dimensionless"/>
<variable initial_value="1" name="theta_27" units="dimensionless"/>
<variable name="K_Concentration" public_interface="out" units="M"/>
<variable name="Na_Concentration" public_interface="out" units="M"/>
<variable name="Cl_Concentration" public_interface="out" units="M"/>
<variable name="Total_Gl_Flux_Apical" public_interface="out" units="umol_per_s"/>
<!--var Concentration_Na: M {pub: out};-->
<!-- Concentration_Na = Na_m+Na_i+Na_s;-->
<!--var K_g: per_second {init: 0.0015, pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>Total_Gl_Flux_Apical</ci>
<apply>
<plus/>
<ci>J_A_GLUT</ci>
<ci>J_Gl_SGLT</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>Na_Concentration</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<ci>theta_15</ci>
</apply>
<ci>J_ENaC</ci>
</apply>
<apply>
<times/>
<ci>theta_3</ci>
<ci>J_NHE3_Na</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>K_Concentration</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_K_NaK</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_NaKCC</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_D_K</ci>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>J_D_Bl_K</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">0</cn>
<ci>J_IK</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Cl_Concentration</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
<apply>
<minus/>
<apply>
<times/>
<apply>
<plus/>
<ci>theta_11</ci>
</apply>
<ci>J_CFTR</ci>
</apply>
<apply>
<times/>
<ci>theta_14</ci>
<ci>J_AE1_Cl</ci>
</apply>
</apply>
</apply>
</apply>
<!--var K_g: per_second {init: 0.0015, pub: out};-->
<apply>
<eq/>
<ci>beta</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2.3</cn>
<ci>HCO3_i</ci>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Na_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<ci>theta_1</ci>
<ci>J_Na_SGLT</ci>
</apply>
<apply>
<times/>
<ci>theta_2</ci>
<ci>J_Na_NaK</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_15</ci>
<ci>J_ENaC</ci>
</apply>
<apply>
<times/>
<ci>theta_16</ci>
<ci>J_NaKCC</ci>
</apply>
<apply>
<times/>
<ci>theta_19</ci>
<ci>J_NBC_Na</ci>
</apply>
<apply>
<times/>
<ci>theta_4</ci>
<ci>J_D_Na</ci>
</apply>
<apply>
<times/>
<ci>theta_5</ci>
<ci>J_D_Bl_Na</ci>
</apply>
<apply>
<times/>
<ci>theta_3</ci>
<ci>J_NHE3_Na</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>glucose_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<minus/>
<apply>
<plus/>
<apply>
<times/>
<ci>theta_26</ci>
<cn cellml:units="dimensionless">1</cn>
<ci>J_A_GLUT</ci>
</apply>
<apply>
<times/>
<ci>theta_6</ci>
<ci>J_Gl_SGLT</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_7</ci>
<cn cellml:units="dimensionless">1</cn>
<ci>J_GLUT</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>K_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>theta_8</ci>
<ci>J_K_NaK</ci>
</apply>
<apply>
<times/>
<ci>theta_17</ci>
<ci>J_NaKCC</ci>
</apply>
<apply>
<times/>
<ci>theta_9</ci>
<ci>J_D_K</ci>
</apply>
<apply>
<times/>
<ci>theta_10</ci>
<ci>J_D_Bl_K</ci>
</apply>
<apply>
<times/>
<ci>theta_27</ci>
<ci>J_IK</ci>
</apply>
<ci>J_BK</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>Cl_i</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>v_cell</ci>
</apply>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<times/>
<apply>
<minus/>
<ci>theta_11</ci>
</apply>
<ci>J_CFTR</ci>
</apply>
<ci>J_CLC2</ci>
</apply>
<apply>
<times/>
<ci>theta_18</ci>
<cn cellml:units="dimensionless">2</cn>
<ci>J_NaKCC</ci>
</apply>
<apply>
<times/>
<ci>theta_12</ci>
<ci>J_D_Cl</ci>
</apply>
<apply>
<times/>
<ci>theta_13</ci>
<ci>J_D_Bl_Cl</ci>
</apply>
<apply>
<times/>
<ci>theta_14</ci>
<ci>J_AE1_Cl</ci>
</apply>
</apply>
</apply>
</apply>
<!--ode(HCO3_i, time) = 1e-9{m3_mole_per_litre_per_umol}/v_cell*theta_20*J_NBC_HCO3+theta_25*J_AE1_HCO3+theta_21*J_buf_C;-->
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>pH_int</ci>
</apply>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="m3_mole_per_litre_per_umol" type="e-notation">1<sep/>-9</cn>
<ci>rho</ci>
</apply>
<apply>
<times/>
<ci>beta</ci>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<ci>theta_20</ci>
<ci>J_NBC_HCO3</ci>
</apply>
<apply>
<times/>
<ci>theta_22</ci>
<ci>J_NHE3_H</ci>
</apply>
</apply>
<apply>
<times/>
<ci>theta_25</ci>
<ci>J_AE1_HCO3</ci>
</apply>
</apply>
</apply>
</apply>
<!--ode(pH_Ext, time) = 1e-9{m3_mole_per_litre_per_umol}*rho/(beta*A_Apical)*(theta_22*J_NHE3_H-theta_25*J_AE1_HCO3-theta_26*J_CFTR_HCO3);-->
<!--ode(pH_s, time) = -1e-9{m3_mole_per_litre_per_umol}*rho/(beta*A_Apical)*theta_20*J_NBC_HCO3;-->
<!--ode(HCO3_m, time) = theta_24*-J_AE1_HCO3+theta_23*J_buf_A;-->
<apply>
<eq/>
<ci>H_int</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pH_int</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>HCO3_i</ci>
<apply>
<times/>
<ci>s</ci>
<ci>P_CO2</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pH_int</ci>
<ci>p_K</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>H_Ext</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pH_Ext</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>H_s</ci>
<apply>
<power/>
<cn cellml:units="dimensionless">10</cn>
<apply>
<minus/>
<ci>pH_s</ci>
</apply>
</apply>
</apply>
</math>
<!-- HCO3_m = s*P_CO2*pow(10{dimensionless}, pH_Ext-p_K);
HCO3_s = s*P_CO2*pow(10{dimensionless}, pH_s-p_K);-->
</component>
<component name="Basolateral_concentrations">
<variable initial_value="0.155" name="Na_s" public_interface="out" units="M"/>
<variable initial_value="0.005" name="glucose_s" public_interface="out" units="M"/>
<variable initial_value="0.119" name="Cl_s" public_interface="out" units="M"/>
<variable initial_value="0.0054" name="K_s" public_interface="out" units="M"/>
<!--var time: second {pub: out};-->
</component>
<component name="Basol_voltage">
<variable initial_value="-0.05" name="v_sc" public_interface="out" units="volt"/>
<variable name="time" public_interface="in" units="second"/>
<variable name="capacitance" public_interface="in" units="uF"/>
<variable name="I_Na_NaK" public_interface="in" units="uA"/>
<variable name="I_K_NaK" public_interface="in" units="uA"/>
<variable name="I_D_P_Na" public_interface="in" units="uA"/>
<variable name="I_D_P_K" public_interface="in" units="uA"/>
<variable name="I_D_P_Cl" public_interface="in" units="uA"/>
<variable name="I_NBC_Na" public_interface="in" units="uA"/>
<variable name="I_NBC_HCO3" public_interface="in" units="uA"/>
<variable name="I_IK" public_interface="in" units="uA"/>
<variable name="I_CLC2" public_interface="in" units="uA"/>
<!--var I_GLUT: uA {pub: in};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<apply>
<diff/>
<bvar>
<ci>time</ci>
</bvar>
<ci>v_sc</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<plus/>
<ci>I_CLC2</ci>
<apply>
<minus/>
<ci>I_Na_NaK</ci>
</apply>
<ci>I_K_NaK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_NBC_Na</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_NBC_HCO3</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_D_P_Na</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_D_P_K</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_D_P_Cl</ci>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="dimensionless">1</cn>
<ci>I_IK</ci>
</apply>
</apply>
</apply>
<ci>capacitance</ci>
</apply>
</apply>
</math>
</component>
<component name="rate_constants">
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="glucose_i" public_interface="in" units="M"/>
<variable name="glucose_m" public_interface="in" units="M"/>
<variable name="k_12" public_interface="out" units="per_second"/>
<variable name="ks_12" public_interface="out" units="per_M2_per_second"/>
<variable name="k_21" public_interface="out" units="per_second"/>
<variable name="k_32" public_interface="out" units="per_second"/>
<variable name="k_25" public_interface="out" units="per_second"/>
<variable name="k_52" public_interface="out" units="per_second"/>
<variable name="k_34" public_interface="out" units="per_second"/>
<variable name="k_43" public_interface="out" units="per_second"/>
<variable name="k_45" public_interface="out" units="per_second"/>
<variable name="k_54" public_interface="out" units="per_second"/>
<variable name="k_56" public_interface="out" units="per_second"/>
<variable name="k_65" public_interface="out" units="per_second"/>
<variable name="k_61" public_interface="out" units="per_second"/>
<variable name="k_16" public_interface="out" units="per_second"/>
<variable name="k0_12" public_interface="in" units="per_M2_per_second"/>
<variable name="k0_21" public_interface="in" units="per_second"/>
<variable name="k0_32" public_interface="in" units="per_second"/>
<variable name="k0_25" public_interface="in" units="per_second"/>
<variable name="k0_52" public_interface="in" units="per_second"/>
<variable name="k0_34" public_interface="in" units="per_second"/>
<variable name="k0_43" public_interface="in" units="per_second"/>
<variable name="k0_45" public_interface="in" units="per_second"/>
<variable name="k0_54" public_interface="in" units="per_M_per_second"/>
<variable name="k0_56" public_interface="in" units="per_second"/>
<variable name="k0_65" public_interface="in" units="per_M2_per_second"/>
<variable name="k0_61" public_interface="in" units="per_second"/>
<variable name="k0_16" public_interface="in" units="per_second"/>
<variable name="delta" public_interface="in" units="dimensionless"/>
<variable name="n" public_interface="in" units="dimensionless"/>
<variable name="z_c" public_interface="in" units="dimensionless"/>
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="alpha_p" public_interface="in" units="dimensionless"/>
<variable name="mu" public_interface="in" units="dimensionless"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>ks_12</ci>
<apply>
<times/>
<ci>k0_12</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>n</ci>
</apply>
<ci>alpha_p</ci>
<ci>mu</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_12</ci>
<apply>
<times/>
<ci>k0_12</ci>
<apply>
<exp/>
<apply>
<times/>
<apply>
<minus/>
<ci>alpha_p</ci>
</apply>
<ci>mu</ci>
</apply>
</apply>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_21</ci>
<apply>
<times/>
<ci>k0_21</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>n</ci>
<ci>z_Na</ci>
<ci>alpha_p</ci>
<ci>mu</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_32</ci>
<ci>k0_32</ci>
</apply>
<apply>
<eq/>
<ci>k_34</ci>
<ci>k0_34</ci>
</apply>
<apply>
<eq/>
<ci>k_43</ci>
<ci>k0_43</ci>
</apply>
<apply>
<eq/>
<ci>k_45</ci>
<ci>k0_45</ci>
</apply>
<apply>
<eq/>
<ci>k_54</ci>
<apply>
<times/>
<ci>k0_54</ci>
<ci>glucose_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>k_25</ci>
<ci>k0_25</ci>
</apply>
<apply>
<eq/>
<ci>k_56</ci>
<ci>k0_56</ci>
</apply>
<apply>
<eq/>
<ci>k_65</ci>
<apply>
<times/>
<ci>k0_65</ci>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_61</ci>
<apply>
<times/>
<ci>k0_61</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_c</ci>
<ci>delta</ci>
<ci>mu</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_16</ci>
<apply>
<times/>
<ci>k0_16</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<ci>z_c</ci>
</apply>
<ci>delta</ci>
<ci>mu</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>k_52</ci>
<ci>k0_52</ci>
</apply>
</math>
<!--k0_54_temp = k0_23*k_34*k_45/(k_43*k_32);-->
</component>
<component name="phenomonological_constants">
<variable name="epsilon" public_interface="out" units="per_second"/>
<variable name="lambda" public_interface="out" units="per_M3_per_second5"/>
<variable name="chi" public_interface="out" units="M"/>
<variable name="alpha" public_interface="out" units="M3"/>
<variable name="beta" public_interface="out" units="M2"/>
<variable name="gamma" public_interface="out" units="M3_per_second"/>
<variable name="phi" public_interface="out" units="M_per_second"/>
<!--var k_12: per_second {pub: in};-->
<variable name="k_21" public_interface="in" units="per_second"/>
<!--var k_23: per_second {pub: in};-->
<variable name="k_32" public_interface="in" units="per_second"/>
<variable name="k_34" public_interface="in" units="per_second"/>
<variable name="k_43" public_interface="in" units="per_second"/>
<variable name="k_45" public_interface="in" units="per_second"/>
<variable name="k_54" public_interface="in" units="per_second"/>
<variable name="k_25" public_interface="in" units="per_second"/>
<variable name="k_52" public_interface="in" units="per_second"/>
<variable name="k_56" public_interface="in" units="per_second"/>
<variable name="k_65" public_interface="in" units="per_second"/>
<variable name="k_61" public_interface="in" units="per_second"/>
<variable name="k_16" public_interface="in" units="per_second"/>
<variable name="ks_12" public_interface="in" units="per_M2_per_second"/>
<variable name="k0_23" public_interface="in" units="per_M_per_second"/>
<!--var C_T: umol {pub: in};-->
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="glucose_m" public_interface="in" units="M"/>
<variable initial_value="40000000" name="n_SGLT" public_interface="out" units="dimensionless"/>
<variable name="J_Na_SGLT" public_interface="out" units="umol_per_s"/>
<variable name="I_Na_SGLT" public_interface="out" units="uA"/>
<variable name="J_Gl_SGLT" public_interface="out" units="umol_per_s"/>
<variable name="J_GL_SGLT1" public_interface="out" units="umol_per_s"/>
<!--var I_Gl_SGLT: uA {pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>lambda</ci>
<apply>
<times/>
<ci>ks_12</ci>
<ci>k0_23</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_43</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>chi</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>lambda</ci>
</apply>
<ci>ks_12</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>lambda</ci>
</apply>
<ci>k0_23</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>lambda</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_54</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_65</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
<ci>k_61</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_56</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
</apply>
<apply>
<times/>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>lambda</ci>
</apply>
<ci>k_16</ci>
<ci>k_21</ci>
<ci>k_65</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_54</ci>
</apply>
<apply>
<times/>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_52</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_45</ci>
<ci>k_52</ci>
</apply>
<apply>
<times/>
<ci>k_32</ci>
<ci>k_43</ci>
<ci>k_52</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>phi</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>lambda</ci>
</apply>
<ci>ks_12</ci>
<ci>k_56</ci>
<ci>k_61</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_34</ci>
<ci>k_45</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_45</ci>
</apply>
<apply>
<times/>
<ci>k_25</ci>
<ci>k_32</ci>
<ci>k_43</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>epsilon</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
<ci>lambda</ci>
</apply>
<ci>ks_12</ci>
<ci>k0_23</ci>
<ci>k_34</ci>
<ci>k_45</ci>
<ci>k_56</ci>
<ci>k_61</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Na_SGLT</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>n_SGLT</ci>
</apply>
<cn cellml:units="per_umol" type="e-notation">6.023<sep/>17</cn>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>epsilon</ci>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>glucose_m</ci>
</apply>
<apply>
<times/>
<ci>phi</ci>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<ci>gamma</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>alpha</ci>
<apply>
<times/>
<ci>beta</ci>
<ci>glucose_m</ci>
</apply>
<apply>
<times/>
<ci>chi</ci>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
<apply>
<times/>
<apply>
<power/>
<ci>Na_m</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>glucose_m</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_Na_SGLT</ci>
<apply>
<times/>
<ci>J_Na_SGLT</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Gl_SGLT</ci>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<cn cellml:units="dimensionless">2</cn>
</apply>
<ci>J_Na_SGLT</ci>
</apply>
</apply>
<!--to change the unit to picomol/hour-->
<apply>
<eq/>
<ci>J_GL_SGLT1</ci>
<apply>
<times/>
<ci>J_Gl_SGLT</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
</math>
<!--I_Gl_SGLT = -J_Gl_SGLT*F;-->
</component>
<component name="CFTR">
<variable name="z_Cl" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="Cl_m" public_interface="in" units="M"/>
<variable name="V_A_Cl" public_interface="out" units="volt"/>
<variable name="V_A_HCO3" public_interface="out" units="volt"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="J_CFTR" public_interface="out" units="umol_per_s"/>
<!-- var J_CFTR1: umol_per_s {pub: out};-->
<variable name="J_CFTR_HCO3" public_interface="out" units="umol_per_s"/>
<variable name="I_CFTR" public_interface="out" units="uA"/>
<variable name="I_CFTR_HCO3" public_interface="out" units="uA"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="1" name="G_CFTR" public_interface="out" units="S_per_m2"/>
<variable name="HCO3_i" public_interface="in" units="M"/>
<variable name="HCO3_m" public_interface="in" units="M"/>
<variable initial_value="-1" name="z_HCO3" public_interface="out" units="dimensionless"/>
<!-- alpha is multipled by area in order to have a dimensionless unit-->
<!--var alpha_CFTR: per_m2 {init: 0.1, pub: out};
J_CFTR1 = J_CFTR*36e8{dimensionless};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_CFTR</ci>
<apply>
<divide/>
<ci>I_CFTR</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Cl_m</ci>
<ci>Cl_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_CFTR</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
</apply>
<ci>G_CFTR</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>V_A_Cl</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_HCO3</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_HCO3</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>HCO3_m</ci>
<ci>HCO3_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_CFTR_HCO3</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">0.25<sep/>6</cn>
<ci>G_CFTR</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>V_A_HCO3</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_CFTR_HCO3</ci>
<apply>
<divide/>
<ci>I_CFTR_HCO3</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="CLC2">
<variable name="z_Cl" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="Cl_s" public_interface="in" units="M"/>
<variable name="V_A_Cl" public_interface="out" units="volt"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="J_CLC2" public_interface="out" units="umol_per_s"/>
<variable name="I_CLC2" public_interface="out" units="uA"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="1" name="G_CLC2" public_interface="out" units="S_per_m2"/>
<!-- alpha is multipled by area in order to have a dimensionless unit-->
<!--var alpha_CFTR: per_m2 {init: 0.1, pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_CLC2</ci>
<apply>
<divide/>
<ci>I_CLC2</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Cl_s</ci>
<ci>Cl_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_CLC2</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
</apply>
<ci>G_CLC2</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_sc</ci>
<ci>V_A_Cl</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="ENaC">
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="V_A_Na" public_interface="out" units="volt"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="J_ENaC" public_interface="out" units="umol_per_s"/>
<variable name="J_ENaC1" public_interface="out" units="umol_per_s"/>
<variable name="I_ENaC" public_interface="out" units="uA"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="2" name="G_ENaC" public_interface="out" units="S_per_m2"/>
<!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
<!--var alpha_ENaC: per_m2 {init: 0.1, pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_ENaC</ci>
<apply>
<divide/>
<ci>I_ENaC</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_ENaC1</ci>
<apply>
<times/>
<ci>J_ENaC</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Na_m</ci>
<ci>Na_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_ENaC</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
</apply>
<ci>G_ENaC</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>V_A_Na</ci>
</apply>
</apply>
</apply>
</math>
</component>
<component name="NaKCC">
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="S_C" public_interface="out" units="M"/>
<variable name="J_NaKCC" public_interface="out" units="umol_per_s"/>
<variable name="I_NaKCC" public_interface="out" units="uA"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="u_NaKCC" public_interface="out" units="per_second"/>
<variable initial_value="157.55" name="gamma_1" public_interface="out" units="per_second2"/>
<variable initial_value="20000000" name="gamma_2" public_interface="out" units="per_M4_per_second2"/>
<variable initial_value="1.03" name="gamma_3" public_interface="out" units="per_second"/>
<variable initial_value="1380000" name="gamma_4" public_interface="out" units="per_M4_per_second"/>
<variable initial_value="1e-9" name="alpha_NaKCC" public_interface="out" units="mole_per_m2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>S_C</ci>
<apply>
<times/>
<ci>Na_i</ci>
<ci>Cl_i</ci>
<ci>Cl_i</ci>
<ci>K_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>u_NaKCC</ci>
<apply>
<divide/>
<apply>
<minus/>
<ci>gamma_1</ci>
<apply>
<times/>
<ci>gamma_2</ci>
<ci>S_C</ci>
</apply>
</apply>
<apply>
<plus/>
<ci>gamma_3</ci>
<apply>
<times/>
<ci>gamma_4</ci>
<ci>S_C</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NaKCC</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
<ci>alpha_NaKCC</ci>
<ci>u_NaKCC</ci>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NaKCC</ci>
<apply>
<times/>
<ci>J_NaKCC</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="A_GLUT2">
<variable initial_value="1000" name="K_12" public_interface="out" units="per_M_per_second"/>
<variable initial_value="12.6" name="K_21" public_interface="out" units="per_second"/>
<variable initial_value="1350" name="K_23" public_interface="out" units="per_second"/>
<variable initial_value="90" name="K_32" public_interface="out" units="per_second"/>
<variable initial_value="12.6" name="K_34" public_interface="out" units="per_second"/>
<variable initial_value="1000" name="K_43" public_interface="out" units="per_M_per_second"/>
<variable initial_value="0.8" name="K_41" public_interface="out" units="per_second"/>
<variable initial_value="12" name="K_14" public_interface="out" units="per_second"/>
<!-- var C_1: M {init: 1e-3, pub: out};
var C_2: M {init: 0, pub: out};-->
<variable name="glucose_i" public_interface="in" units="M"/>
<variable name="glucose_m" public_interface="in" units="M"/>
<!-- var CS_1: M {pub: out};
var CS_2: M {pub: out};-->
<variable name="K_S" public_interface="out" units="M"/>
<variable name="J_A_GLUT" public_interface="out" units="umol_per_s"/>
<variable name="C_tot" public_interface="out" units="umol"/>
<variable initial_value="420000000" name="n_GLUT" public_interface="out" units="dimensionless"/>
<variable name="D" public_interface="out" units="per_second"/>
<variable name="J_A_GLUT1" public_interface="out" units="umol_per_s"/>
<!--var time: second {init: 0, pub: out};-->
<variable name="v_cell" public_interface="in" units="m3"/>
<!--var I_GLUT: uA {pub: out};-->
<!--var n_GLUT: dimensionless {init: 13.2e7};
var F: C_per_mol {pub: in};-->
<!-- D = sel
case (glucose_m > 0.02{M}):
(1{dimensionless}+glucose_i/K_S)*(K_32*glucose_m/K_S+K_41)+(1{dimensionless}+glucose_m/K_S)*(K_23*glucose_i/K_S+K_14);
case glucose_m < 0.02{M}:
1e100{per_second};
endsel;-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>C_tot</ci>
<apply>
<divide/>
<ci>n_GLUT</ci>
<cn cellml:units="per_umole" type="e-notation">6.022<sep/>17</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>K_S</ci>
<apply>
<divide/>
<ci>K_21</ci>
<ci>K_12</ci>
</apply>
</apply>
<!-- CS_1 = C_tot* ((K_32*(S_2)/K_S) + K_41)*(S_1/K_S);
CS_2 = C_tot* ((K_23*(S_1)/K_S) + K_14)*(S_2/K_S);-->
<apply>
<eq/>
<ci>J_A_GLUT</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>C_tot</ci>
<ci>D</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_32</ci>
<ci>K_14</ci>
<ci>glucose_m</ci>
</apply>
<ci>K_S</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>K_23</ci>
<ci>K_41</ci>
<ci>glucose_i</ci>
</apply>
<ci>K_S</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_A_GLUT1</ci>
<apply>
<times/>
<ci>J_A_GLUT</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>D</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>glucose_m</ci>
<ci>K_S</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_23</ci>
<ci>glucose_i</ci>
</apply>
<ci>K_S</ci>
</apply>
<ci>K_14</ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>glucose_i</ci>
<ci>K_S</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_32</ci>
<ci>glucose_m</ci>
</apply>
<ci>K_S</ci>
</apply>
<ci>K_41</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="GLUT2">
<variable initial_value="1000" name="K_12" public_interface="out" units="per_M_per_second"/>
<variable initial_value="12.6" name="K_21" public_interface="out" units="per_second"/>
<variable initial_value="1350" name="K_23" public_interface="out" units="per_second"/>
<variable initial_value="90" name="K_32" public_interface="out" units="per_second"/>
<variable initial_value="12.6" name="K_34" public_interface="out" units="per_second"/>
<variable initial_value="1000" name="K_43" public_interface="out" units="per_M_per_second"/>
<variable initial_value="0.8" name="K_41" public_interface="out" units="per_second"/>
<variable initial_value="12" name="K_14" public_interface="out" units="per_second"/>
<!-- var C_1: M {init: 1e-3, pub: out};
var C_2: M {init: 0, pub: out};-->
<variable name="glucose_i" public_interface="in" units="M"/>
<variable name="glucose_s" public_interface="in" units="M"/>
<!-- var CS_1: M {pub: out};
var CS_2: M {pub: out};-->
<variable name="K_S" public_interface="out" units="M"/>
<variable name="J_GLUT" public_interface="out" units="umol_per_s"/>
<variable name="C_tot" public_interface="out" units="umol"/>
<variable initial_value="140000000" name="n_GLUT" public_interface="out" units="dimensionless"/>
<variable name="D" public_interface="out" units="per_second"/>
<variable name="J_GLUT1" public_interface="out" units="umol_per_s"/>
<!--var time: second {init: 0, pub: out};-->
<variable name="v_cell" public_interface="in" units="m3"/>
<!--var I_GLUT: uA {pub: out};-->
<!--var n_GLUT: dimensionless {init: 13.2e7};
var F: C_per_mol {pub: in};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>C_tot</ci>
<apply>
<divide/>
<ci>n_GLUT</ci>
<cn cellml:units="per_umole" type="e-notation">6.022<sep/>17</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>K_S</ci>
<apply>
<divide/>
<ci>K_21</ci>
<ci>K_12</ci>
</apply>
</apply>
<!-- CS_1 = C_tot* ((K_32*(S_2)/K_S) + K_41)*(S_1/K_S);
CS_2 = C_tot* ((K_23*(S_1)/K_S) + K_14)*(S_2/K_S);-->
<apply>
<eq/>
<ci>J_GLUT</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>C_tot</ci>
<ci>D</ci>
</apply>
<apply>
<minus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_23</ci>
<ci>K_41</ci>
<ci>glucose_i</ci>
</apply>
<ci>K_S</ci>
</apply>
<apply>
<divide/>
<apply>
<times/>
<ci>K_32</ci>
<ci>K_14</ci>
<ci>glucose_s</ci>
</apply>
<ci>K_S</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_GLUT1</ci>
<apply>
<times/>
<ci>J_GLUT</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>D</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>glucose_i</ci>
<ci>K_S</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_32</ci>
<ci>glucose_s</ci>
</apply>
<ci>K_S</ci>
</apply>
<ci>K_41</ci>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>glucose_s</ci>
<ci>K_S</ci>
</apply>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<apply>
<times/>
<ci>K_23</ci>
<ci>glucose_i</ci>
</apply>
<ci>K_S</ci>
</apply>
<ci>K_14</ci>
</apply>
</apply>
</apply>
</apply>
</math>
</component>
<component name="NaK_channel">
<variable name="time" public_interface="in" units="second"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="K_s" public_interface="in" units="M"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="TO_NaK" public_interface="out" units="per_second"/>
<variable initial_value="5000000" name="n_NaK" public_interface="out" units="dimensionless"/>
<variable initial_value="0.01" name="K_Na_i_half" public_interface="out" units="M"/>
<variable initial_value="0.0015" name="K_K_s_half" public_interface="out" units="M"/>
<variable name="I_Na_NaK" public_interface="out" units="uA"/>
<variable name="I_K_NaK" public_interface="out" units="uA"/>
<variable name="J_NaK" public_interface="out" units="umol_per_s"/>
<variable name="J_Na_NaK" public_interface="out" units="umol_per_s"/>
<variable name="J_K_NaK" public_interface="out" units="umol_per_s"/>
<!--var J_K_NaK: umol_per_s {pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>I_K_NaK</ci>
<apply>
<times/>
<ci>J_K_NaK</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_Na_NaK</ci>
<apply>
<times/>
<ci>J_Na_NaK</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NaK</ci>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>TO_NaK</ci>
<ci>n_NaK</ci>
</apply>
<cn cellml:units="per_umol" type="e-notation">6.023<sep/>17</cn>
</apply>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">1.36</cn>
</apply>
<cn cellml:units="dimensionless">1.13</cn>
</apply>
<apply>
<plus/>
<apply>
<power/>
<ci>K_Na_i_half</ci>
<cn cellml:units="dimensionless">1.36</cn>
</apply>
<apply>
<power/>
<ci>Na_i</ci>
<cn cellml:units="dimensionless">1.36</cn>
</apply>
</apply>
</apply>
<cn cellml:units="dimensionless">1.3</cn>
<ci>K_s</ci>
</apply>
<apply>
<plus/>
<ci>K_K_s_half</ci>
<ci>K_s</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>TO_NaK</ci>
<apply>
<plus/>
<apply>
<minus/>
<apply>
<minus/>
<apply>
<times/>
<cn cellml:units="per_volt4_per_second" type="e-notation">5.46<sep/>+4</cn>
<apply>
<power/>
<ci>v_sc</ci>
<cn cellml:units="dimensionless">4</cn>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_volt3_per_second" type="e-notation">2.43<sep/>+2</cn>
<apply>
<power/>
<ci>v_sc</ci>
<cn cellml:units="dimensionless">3</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_volt2_per_second" type="e-notation">2<sep/>+3</cn>
<apply>
<power/>
<ci>v_sc</ci>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
<apply>
<times/>
<cn cellml:units="per_volt_per_second">160</cn>
<ci>v_sc</ci>
</apply>
<cn cellml:units="per_second">51</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Na_NaK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">3</cn>
<ci>J_NaK</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_K_NaK</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">2</cn>
<ci>J_NaK</ci>
</apply>
</apply>
</math>
</component>
<component name="NBC">
<variable name="time" public_interface="in" units="second"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="v_cell" public_interface="in" units="m3"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="Na_s" public_interface="in" units="M"/>
<variable name="HCO3_i" public_interface="in" units="M"/>
<variable name="HCO3_s" public_interface="in" units="M"/>
<variable initial_value="0.08" name="g_NBC" public_interface="out" units="S_per_m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<!--kelvin {init: 300};-->
<variable initial_value="2" name="n" units="dimensionless"/>
<!--var eps: M {init: 1e-12};-->
<variable name="E_NBC" public_interface="out" units="volt"/>
<variable name="J_NBC_Na" public_interface="out" units="umol_per_s"/>
<variable name="J_NBC_HCO3" public_interface="out" units="umol_per_s"/>
<variable name="I_NBC_Na" public_interface="out" units="uA"/>
<variable name="I_NBC_HCO3" public_interface="out" units="uA"/>
<variable initial_value="1e-18" name="eps" units="M"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>E_NBC</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<ci>n</ci>
</apply>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<times/>
<apply>
<divide/>
<ci>Na_s</ci>
<ci>Na_i</ci>
</apply>
<apply>
<power/>
<apply>
<divide/>
<ci>HCO3_s</ci>
<ci>HCO3_i</ci>
</apply>
<cn cellml:units="dimensionless">2</cn>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NBC_Na</ci>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
<ci>g_NBC</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_sc</ci>
<ci>E_NBC</ci>
</apply>
</apply>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NBC_HCO3</ci>
<apply>
<times/>
<cn cellml:units="dimensionless">3.0</cn>
<ci>J_NBC_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NBC_Na</ci>
<apply>
<times/>
<ci>J_NBC_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NBC_HCO3</ci>
<apply>
<times/>
<ci>J_NBC_HCO3</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="NHE3">
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="v_cell" public_interface="in" units="m3"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="H_Ext" public_interface="in" units="M"/>
<variable name="H_int" public_interface="in" units="M"/>
<variable initial_value="0" name="NH4_ext" public_interface="out" units="M"/>
<variable initial_value="0" name="NH4_int" public_interface="out" units="M"/>
<variable name="J_NHE3_Na" public_interface="out" units="umol_per_s"/>
<variable name="J_NHE3_H1" public_interface="out" units="umol_per_s"/>
<variable name="J_NHE3_H" public_interface="out" units="umol_per_s"/>
<variable name="J_NHE3_NH4" public_interface="out" units="umol_per_s"/>
<variable name="J_NHE3_Na_Max" public_interface="out" units="umol_per_s"/>
<variable name="I_NHE3_H" public_interface="out" units="uA"/>
<variable name="I_NHE3_Na" public_interface="out" units="uA"/>
<variable initial_value="0.001" name="x_T" units="M"/>
<variable name="sigma" units="per_second"/>
<variable name="P_Na" units="per_second"/>
<variable name="P_H" units="per_second"/>
<variable name="P_NH4" units="per_second"/>
<variable initial_value="1.6" name="P0_Na" units="per_second"/>
<variable initial_value="0.48" name="P0_H" units="per_second"/>
<variable initial_value="1.6" name="P0_NH4" units="per_second"/>
<variable initial_value="0.03" name="K_Na" units="M"/>
<variable initial_value="7.2e-8" name="K_H" units="M"/>
<variable initial_value="0.027" name="K_NH4" units="M"/>
<variable initial_value="1e-9" name="K_I" units="M"/>
<variable initial_value="0" name="f_m" units="dimensionless"/>
<variable initial_value="2" name="f_M" units="dimensionless"/>
<variable name="alpha_ext_Na" units="dimensionless"/>
<variable name="alpha_int_Na" units="dimensionless"/>
<variable name="beta_ext_H" units="dimensionless"/>
<variable name="beta_int_H" units="dimensionless"/>
<variable name="gamma_ext_NH4" units="dimensionless"/>
<variable name="gamma_int_NH4" units="dimensionless"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<!-- var pH_Ext: dimensionless;
var pH_int: dimensionless;
var p_K: dimensionless {init: 6.1, pub: out};
var P_CO2: dimensionless {init: 40, pub: out};
var s: dimensionless {init: 0.03e-3, pub: out};
var HCO3_i: M {pub: in};
var HCO3_m: M {pub: in};-->
<!-- pH_Ext = p_K+log(HCO3_m/(s*P_CO2));
pH_int = p_K+log(HCO3_i/(s*P_CO2));-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>alpha_ext_Na</ci>
<apply>
<divide/>
<ci>Na_m</ci>
<ci>K_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>alpha_int_Na</ci>
<apply>
<divide/>
<ci>Na_i</ci>
<ci>K_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_ext_H</ci>
<apply>
<divide/>
<ci>H_Ext</ci>
<ci>K_H</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_int_H</ci>
<apply>
<divide/>
<ci>H_int</ci>
<ci>K_H</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma_ext_NH4</ci>
<apply>
<divide/>
<ci>NH4_ext</ci>
<ci>K_NH4</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma_int_NH4</ci>
<apply>
<divide/>
<ci>NH4_int</ci>
<ci>K_NH4</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>P_Na</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>P0_Na</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>f_M</ci>
<ci>H_int</ci>
</apply>
<apply>
<times/>
<ci>f_m</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci>H_int</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_H</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>P0_H</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>f_M</ci>
<ci>H_int</ci>
</apply>
<apply>
<times/>
<ci>f_m</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci>H_int</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>P_NH4</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>P0_NH4</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>f_M</ci>
<ci>H_int</ci>
</apply>
<apply>
<times/>
<ci>f_m</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<plus/>
<ci>H_int</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>sigma</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>alpha_ext_Na</ci>
<ci>beta_ext_H</ci>
<ci>gamma_ext_NH4</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>alpha_int_Na</ci>
</apply>
<apply>
<times/>
<ci>P_H</ci>
<ci>beta_int_H</ci>
</apply>
<apply>
<times/>
<ci>P_NH4</ci>
<ci>gamma_int_NH4</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>alpha_int_Na</ci>
<ci>beta_int_H</ci>
<ci>gamma_int_NH4</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>alpha_ext_Na</ci>
</apply>
<apply>
<times/>
<ci>P_H</ci>
<ci>beta_ext_H</ci>
</apply>
<apply>
<times/>
<ci>P_NH4</ci>
<ci>gamma_ext_NH4</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
</apply>
<ci>v_cell</ci>
<ci>x_T</ci>
</apply>
<ci>sigma</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>P_H</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_int_Na</ci>
<ci>beta_ext_H</ci>
</apply>
<apply>
<times/>
<ci>alpha_ext_Na</ci>
<ci>beta_int_H</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>P_Na</ci>
<ci>P_NH4</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_int_Na</ci>
<ci>gamma_ext_NH4</ci>
</apply>
<apply>
<times/>
<ci>alpha_ext_Na</ci>
<ci>gamma_int_NH4</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_H</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
</apply>
<ci>v_cell</ci>
<ci>x_T</ci>
</apply>
<ci>sigma</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>P_H</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_ext_Na</ci>
<ci>beta_int_H</ci>
</apply>
<apply>
<times/>
<ci>alpha_int_Na</ci>
<ci>beta_ext_H</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>P_H</ci>
<ci>P_NH4</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>beta_int_H</ci>
<ci>gamma_ext_NH4</ci>
</apply>
<apply>
<times/>
<ci>beta_ext_H</ci>
<ci>gamma_int_NH4</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_H1</ci>
<apply>
<times/>
<ci>J_NHE3_H</ci>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_NH4</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>x_T</ci>
<ci>sigma</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_Na</ci>
<ci>P_NH4</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>alpha_ext_Na</ci>
<ci>gamma_int_NH4</ci>
</apply>
<apply>
<times/>
<ci>alpha_int_Na</ci>
<ci>gamma_ext_NH4</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<ci>P_H</ci>
<ci>P_NH4</ci>
<apply>
<minus/>
<apply>
<times/>
<ci>beta_ext_H</ci>
<ci>gamma_int_NH4</ci>
</apply>
<apply>
<times/>
<ci>gamma_ext_NH4</ci>
<ci>beta_int_H</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_NHE3_Na_Max</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>x_T</ci>
<ci>P_Na</ci>
<ci>P_H</ci>
</apply>
<apply>
<plus/>
<ci>P_Na</ci>
<ci>P_H</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NHE3_H</ci>
<apply>
<times/>
<ci>J_NHE3_H</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_NHE3_Na</ci>
<apply>
<times/>
<ci>J_NHE3_Na</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="AE1">
<!--var A_Apical: m2 {pub: in};-->
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="v_cell" public_interface="in" units="m3"/>
<variable name="HCO3_i" public_interface="in" units="M"/>
<variable name="HCO3_m" public_interface="in" units="M"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="Cl_m" public_interface="in" units="M"/>
<variable name="J_AE1_HCO3" public_interface="out" units="umol_per_s"/>
<!--var J_AE1_HCO31: umol_per_s {pub: out};-->
<variable name="J_AE1_Cl" public_interface="out" units="umol_per_s"/>
<!--var J_AE1_Cl1: umol_per_s {pub: out};-->
<variable name="I_AE1_HCO3" public_interface="out" units="uA"/>
<variable name="I_AE1_Cl" public_interface="out" units="uA"/>
<variable name="J_HCO3_influx" public_interface="out" units="umol_per_s"/>
<variable name="J_Cl_influx" public_interface="out" units="umol_per_s"/>
<variable initial_value="0.198" name="K_HCO3_m" units="M"/>
<variable initial_value="0.198" name="K_HCO3_i" units="M"/>
<variable initial_value="0.05" name="K_Cl_s" units="M"/>
<variable initial_value="0.05" name="K_Cl_i" units="M"/>
<variable initial_value="1247" name="P_HCO3_m" units="per_second"/>
<variable initial_value="135" name="P_HCO3_i" units="per_second"/>
<variable initial_value="562" name="P_Cl_s" units="per_second"/>
<variable initial_value="61" name="P_Cl_i" units="per_second"/>
<variable name="beta_ext" units="dimensionless"/>
<variable name="beta_int" units="dimensionless"/>
<variable name="gamma_ext" units="dimensionless"/>
<variable name="gamma_int" units="dimensionless"/>
<variable name="sigma" units="per_second"/>
<variable initial_value="0.001" name="x_Tmax" units="M"/>
<variable initial_value="0.172" name="K_I" units="M"/>
<variable name="x_T" units="M"/>
<variable name="x_ext" units="M"/>
<variable name="x_int" units="M"/>
<variable name="Jo_bm" units="umol_per_s"/>
<variable name="Ji_bm" units="umol_per_s"/>
<variable name="Js_bm" units="umol_per_s"/>
<variable name="Jo_cm" units="umol_per_s"/>
<variable name="Ji_cm" units="umol_per_s"/>
<variable name="Js_cm" units="umol_per_s"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>x_T</ci>
<apply>
<divide/>
<ci>x_Tmax</ci>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<divide/>
<ci>HCO3_i</ci>
<ci>K_I</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_ext</ci>
<apply>
<divide/>
<ci>HCO3_m</ci>
<ci>K_HCO3_m</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>beta_int</ci>
<apply>
<divide/>
<ci>HCO3_i</ci>
<ci>K_HCO3_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma_ext</ci>
<apply>
<divide/>
<ci>Cl_m</ci>
<ci>K_Cl_s</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>gamma_int</ci>
<apply>
<divide/>
<ci>Cl_i</ci>
<ci>K_Cl_i</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>sigma</ci>
<apply>
<plus/>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>beta_ext</ci>
<ci>gamma_ext</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
<apply>
<times/>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<ci>beta_int</ci>
<ci>gamma_int</ci>
</apply>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_m</ci>
<ci>beta_ext</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_s</ci>
<ci>gamma_ext</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_AE1_HCO3</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>9</cn>
<ci>v_cell</ci>
<ci>x_Tmax</ci>
</apply>
<ci>sigma</ci>
</apply>
<apply>
<minus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
<ci>P_Cl_s</ci>
<ci>gamma_ext</ci>
</apply>
<apply>
<times/>
<ci>P_HCO3_m</ci>
<ci>beta_ext</ci>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
</apply>
<!--J_AE1_HCO31 = J_AE1_HCO3*36e8{dimensionless};
J_AE1_Cl1 = -J_AE1_HCO31;-->
<apply>
<eq/>
<ci>J_AE1_Cl</ci>
<ci>J_AE1_HCO3</ci>
</apply>
<apply>
<eq/>
<ci>I_AE1_HCO3</ci>
<apply>
<times/>
<ci>J_AE1_HCO3</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_AE1_Cl</ci>
<apply>
<times/>
<ci>J_AE1_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>x_ext</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>x_T</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
<ci>sigma</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>x_int</ci>
<apply>
<divide/>
<apply>
<times/>
<ci>x_T</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_m</ci>
<ci>beta_ext</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_s</ci>
<ci>gamma_ext</ci>
</apply>
</apply>
</apply>
<ci>sigma</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_HCO3_influx</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>x_T</ci>
<ci>sigma</ci>
</apply>
<ci>P_HCO3_m</ci>
<ci>beta_ext</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_Cl_influx</ci>
<apply>
<times/>
<apply>
<divide/>
<ci>x_T</ci>
<ci>sigma</ci>
</apply>
<ci>P_Cl_s</ci>
<ci>gamma_ext</ci>
<apply>
<plus/>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>beta_int</ci>
</apply>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>gamma_int</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Jo_bm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_m</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_i</ci>
</apply>
<apply>
<divide/>
<ci>K_HCO3_i</ci>
<apply>
<times/>
<ci>P_HCO3_i</ci>
<ci>HCO3_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ji_bm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_m</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_i</ci>
</apply>
<apply>
<divide/>
<ci>K_HCO3_m</ci>
<apply>
<times/>
<ci>P_HCO3_m</ci>
<ci>HCO3_m</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Js_bm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_m</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_HCO3_i</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Jo_cm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_s</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_i</ci>
</apply>
<apply>
<divide/>
<ci>K_Cl_i</ci>
<apply>
<times/>
<ci>P_Cl_i</ci>
<ci>Cl_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Ji_cm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_s</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_i</ci>
</apply>
<apply>
<divide/>
<ci>K_Cl_s</ci>
<apply>
<times/>
<ci>P_Cl_s</ci>
<ci>Cl_m</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>Js_cm</ci>
<apply>
<power/>
<apply>
<times/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>x_T</ci>
</apply>
<apply>
<plus/>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_s</ci>
</apply>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<ci>P_Cl_i</ci>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
</apply>
</apply>
</apply>
</math>
</component>
<!-- def comp Hydration as
var time: second {pub: in};
var v_cell: m3 {pub: in};
var CO2_i: M {init: 1e-3, pub: out};
var CO2_s: M {init: 1.2e-3, pub: out};
var CO2_m: M {init: 1.2e-3, pub: out};
var H_int: M {pub: in};
var H_Ext: M {pub: in};
var HCO3_i: M {pub: in};
var HCO3_m: M {pub: in};
var p_CO2: per_second {init: 5e-6, pub: out};
var k_f: per_second {init: 0.03, pub: out};
var k_r: per_M_per_second {init: 20e-3, pub: out};
var J_CDF_A: umol_per_s {pub: out};
var J_CDF_B: umol_per_s {pub: out};
var J_buf_C: umol_per_s {pub: out};
var J_buf_C1: umol_per_s {pub: out};
var J_buf_A: umol_per_s {pub: out};
var J_D_P_CO2: umol_per_s {pub: in};-->
<!-- J_CDF_A = -1e9{dimensionless}*v_cell*p_CO2*(CO2_i-CO2_m);
J_CDF_B = 1e9{dimensionless}*v_cell*p_CO2*(CO2_i-CO2_s);
J_buf_C = 1e9{dimensionless}*v_cell*(k_f*CO2_i-k_r*HCO3_i*H_int);
J_buf_C1 = J_buf_C*36e8{dimensionless};
J_buf_A = 1e9{dimensionless}*v_cell*(k_f*CO2_m-k_r*HCO3_m*H_Ext);
ode(CO2_i, time) = -J_buf_C+J_CDF_A+J_CDF_B;
ode(CO2_m, time) = -J_buf_A-J_CDF_A-0{dimensionless}*J_D_P_CO2;
enddef;-->
<component name="IK">
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="K_s" public_interface="in" units="M"/>
<variable name="Ca_i" public_interface="in" units="M"/>
<variable name="V_A_IK" public_interface="out" units="volt"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="J_IK" public_interface="out" units="umol_per_s"/>
<variable name="I_IK" public_interface="out" units="uA"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="1.7" name="G_IK" public_interface="out" units="S_per_m2"/>
<!--var P_IK: dimensionless {pub: out};-->
<variable initial_value="0.0003" name="K_s_IK" public_interface="out" units="M"/>
<variable initial_value="1.7" name="eta" public_interface="out" units="dimensionless"/>
<!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
<!--var alpha_IK: per_m2 {init: 1, pub: out};-->
<!--P_IK = 1{dimensionless}/(1{dimensionless}+pow(K_s_IK/Ca_i, eta));-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>V_A_IK</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>K_s</ci>
<ci>K_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_IK</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
</apply>
<ci>G_IK</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_sc</ci>
<ci>V_A_IK</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_IK</ci>
<apply>
<divide/>
<ci>I_IK</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="BK">
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="K_m" public_interface="in" units="M"/>
<variable name="Ca_i" public_interface="in" units="M"/>
<variable name="V_A_BK" public_interface="out" units="volt"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="J_BK" public_interface="out" units="umol_per_s"/>
<variable name="I_BK" public_interface="out" units="uA"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable initial_value="1" name="G_BK" public_interface="out" units="S_per_m2"/>
<variable name="P_BK" public_interface="out" units="dimensionless"/>
<variable initial_value="0.0003" name="K_s_BK" public_interface="out" units="M"/>
<variable initial_value="1.7" name="eta" public_interface="out" units="dimensionless"/>
<!-- alpha is not multipled by area in order to have a dimensionless unit, so the unit is actually 1/m2, it's necessary to the same for CFTR as well-->
<!--var alpha_BK: per_m2 {init: 0.1, pub: out};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>P_BK</ci>
<apply>
<divide/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<plus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<power/>
<apply>
<divide/>
<ci>K_s_BK</ci>
<ci>Ca_i</ci>
</apply>
<ci>eta</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>V_A_BK</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>K_m</ci>
<ci>K_i</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_BK</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>+6</cn>
</apply>
<ci>G_BK</ci>
<ci>P_BK</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>V_A_BK</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_BK</ci>
<apply>
<divide/>
<ci>I_BK</ci>
<ci>F</ci>
</apply>
</apply>
</math>
</component>
<component name="Apical_Diffusion">
<variable name="v_mc" public_interface="in" units="volt"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="K_m" public_interface="in" units="M"/>
<variable name="Cl_m" public_interface="in" units="M"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="z_Cl" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="P_Na" public_interface="in" units="m_per_s"/>
<variable name="P_k" public_interface="in" units="m_per_s"/>
<variable name="P_Cl" public_interface="in" units="m_per_s"/>
<variable name="J_D_Na" public_interface="out" units="umol_per_s"/>
<variable name="J_D_NA" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Na" public_interface="out" units="uA"/>
<variable name="J_D_K" public_interface="out" units="umol_per_s"/>
<variable name="I_D_K" public_interface="out" units="uA"/>
<variable name="J_D_Cl" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Cl" public_interface="out" units="uA"/>
<variable name="J_D_Apical" public_interface="out" units="umol_per_s"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_D_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Na</ci>
<ci>F</ci>
<ci>z_Na</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Na_m</ci>
<apply>
<times/>
<ci>Na_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_NA</ci>
<apply>
<times/>
<cn cellml:units="dimensionless" type="e-notation">36<sep/>8</cn>
<ci>J_D_Na</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Na</ci>
<apply>
<times/>
<ci>J_D_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_k</ci>
<ci>F</ci>
<ci>z_k</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>K_m</ci>
<apply>
<times/>
<ci>K_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_K</ci>
<apply>
<times/>
<ci>J_D_K</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Cl</ci>
<ci>F</ci>
<ci>z_Cl</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Cl_m</ci>
<apply>
<times/>
<ci>Cl_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>v_mc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Apical</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Cl</ci>
<apply>
<times/>
<ci>J_D_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Apical</ci>
<apply>
<plus/>
<ci>J_D_Na</ci>
<ci>J_D_K</ci>
<ci>J_D_Cl</ci>
</apply>
</apply>
</math>
</component>
<component name="Basolateral_Diffusion">
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="Na_s" public_interface="in" units="M"/>
<variable name="K_s" public_interface="in" units="M"/>
<variable name="Cl_s" public_interface="in" units="M"/>
<variable name="Na_i" public_interface="in" units="M"/>
<variable name="K_i" public_interface="in" units="M"/>
<variable name="Cl_i" public_interface="in" units="M"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="z_Cl" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="A_Basol" public_interface="in" units="m2"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="P_Bl_Na" public_interface="in" units="m_per_s"/>
<variable name="P_Bl_k" public_interface="in" units="m_per_s"/>
<variable name="P_Bl_Cl" public_interface="in" units="m_per_s"/>
<variable name="J_D_Bl_Na" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Bl_Na" public_interface="out" units="uA"/>
<variable name="J_D_Bl_K" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Bl_K" public_interface="out" units="uA"/>
<variable name="J_D_Bl_Cl" public_interface="out" units="umol_per_s"/>
<variable name="I_D_Bl_Cl" public_interface="out" units="uA"/>
<variable name="J_D_Basolateral" public_interface="out" units="umol_per_s"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>J_D_Bl_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Bl_Na</ci>
<ci>F</ci>
<ci>z_Na</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Na_s</ci>
<apply>
<times/>
<ci>Na_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Basol</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Bl_Na</ci>
<apply>
<times/>
<ci>J_D_Bl_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Bl_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Bl_k</ci>
<ci>F</ci>
<ci>z_k</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>K_s</ci>
<apply>
<times/>
<ci>K_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Basol</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Bl_K</ci>
<apply>
<times/>
<ci>J_D_Bl_K</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Bl_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="litre_umol_per_mol_per_m3" type="e-notation">1<sep/>9</cn>
</apply>
<ci>P_Bl_Cl</ci>
<ci>F</ci>
<ci>z_Cl</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
<apply>
<minus/>
<ci>Cl_s</ci>
<apply>
<times/>
<ci>Cl_i</ci>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
</apply>
<apply>
<minus/>
<cn cellml:units="dimensionless">1</cn>
<apply>
<exp/>
<apply>
<divide/>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
<ci>v_sc</ci>
</apply>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
</apply>
</apply>
</apply>
</apply>
<ci>A_Basol</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_Bl_Cl</ci>
<apply>
<times/>
<ci>J_D_Bl_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Basolateral</ci>
<apply>
<plus/>
<ci>J_D_Bl_Na</ci>
<ci>J_D_Bl_K</ci>
<ci>J_D_Bl_Cl</ci>
</apply>
</apply>
</math>
</component>
<component name="Paracellular_voltage">
<variable name="v_ms" public_interface="out" units="volt"/>
<variable name="v_sc" public_interface="in" units="volt"/>
<variable name="v_mc" public_interface="in" units="volt"/>
<!-- var time: second {pub: in};
var capacitance: uF {pub: in};
var I_D_P_Na: uA {pub: in};
var I_D_P_K: uA {pub: in};
var I_D_P_Cl: uA {pub: in};-->
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>v_ms</ci>
<apply>
<minus/>
<ci>v_mc</ci>
<ci>v_sc</ci>
</apply>
</apply>
</math>
<!--ode(v_ms, time) = (I_D_P_Na+I_D_P_K+I_D_P_Cl)/capacitance;-->
</component>
<component name="Paracellular_Diffusion">
<variable name="v_ms" public_interface="in" units="volt"/>
<variable name="Na_m" public_interface="in" units="M"/>
<variable name="K_m" public_interface="in" units="M"/>
<variable name="Cl_m" public_interface="in" units="M"/>
<variable name="Na_s" public_interface="in" units="M"/>
<variable name="K_s" public_interface="in" units="M"/>
<variable name="Cl_s" public_interface="in" units="M"/>
<variable name="F" public_interface="in" units="C_per_mol"/>
<variable name="z_Na" public_interface="in" units="dimensionless"/>
<variable name="z_k" public_interface="in" units="dimensionless"/>
<variable name="z_Cl" public_interface="in" units="dimensionless"/>
<variable name="R" public_interface="in" units="J_per_K_per_mol"/>
<variable name="A_Apical" public_interface="in" units="m2"/>
<variable name="T" public_interface="in" units="kelvin"/>
<variable name="P_P_Na" public_interface="in" units="m_per_s"/>
<variable name="P_P_k" public_interface="in" units="m_per_s"/>
<variable name="P_P_Cl" public_interface="in" units="m_per_s"/>
<variable name="J_D_P_Na" public_interface="out" units="umol_per_s"/>
<variable name="I_D_P_Na" public_interface="out" units="uA"/>
<variable name="J_D_P_K" public_interface="out" units="umol_per_s"/>
<variable name="I_D_P_K" public_interface="out" units="uA"/>
<variable name="J_D_P_Cl" public_interface="out" units="umol_per_s"/>
<variable name="I_D_P_Cl" public_interface="out" units="uA"/>
<variable name="J_D_Paracellular" public_interface="out" units="umol_per_s"/>
<variable name="V_P_Na" public_interface="out" units="volt"/>
<variable name="V_P_K" public_interface="out" units="volt"/>
<variable name="V_P_Cl" public_interface="out" units="volt"/>
<variable initial_value="5" name="G_P_Na" public_interface="out" units="S_per_m2"/>
<variable initial_value="3" name="G_P_K" public_interface="out" units="S_per_m2"/>
<variable initial_value="5" name="G_P_Cl" public_interface="out" units="S_per_m2"/>
<math xmlns="http://www.w3.org/1998/Math/MathML">
<apply>
<eq/>
<ci>V_P_Na</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Na_m</ci>
<ci>Na_s</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_P_Na</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
</apply>
<ci>G_P_Na</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_ms</ci>
<ci>V_P_Na</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_P_Na</ci>
<apply>
<divide/>
<ci>I_D_P_Na</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_P_K</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_k</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>K_m</ci>
<ci>K_s</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_P_K</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
</apply>
<ci>G_P_K</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_ms</ci>
<ci>V_P_K</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_P_K</ci>
<apply>
<divide/>
<ci>I_D_P_K</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>V_P_Cl</ci>
<apply>
<times/>
<apply>
<divide/>
<apply>
<times/>
<ci>R</ci>
<ci>T</ci>
</apply>
<apply>
<times/>
<ci>z_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<ln/>
<apply>
<divide/>
<ci>Cl_m</ci>
<ci>Cl_s</ci>
</apply>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>I_D_P_Cl</ci>
<apply>
<times/>
<apply>
<minus/>
<cn cellml:units="dimensionless" type="e-notation">1<sep/>6</cn>
</apply>
<ci>G_P_Cl</ci>
<ci>A_Apical</ci>
<apply>
<minus/>
<ci>v_ms</ci>
<ci>V_P_Cl</ci>
</apply>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_P_Cl</ci>
<apply>
<divide/>
<ci>I_D_P_Cl</ci>
<ci>F</ci>
</apply>
</apply>
<apply>
<eq/>
<ci>J_D_Paracellular</ci>
<apply>
<plus/>
<ci>J_D_P_Na</ci>
<ci>J_D_P_K</ci>
<ci>J_D_P_Cl</ci>
</apply>
</apply>
</math>
</component>
<!-- def comp interstitium_concentration as
var glucose_s: M {init: 0.005, pub: out};-->
<!-- //var J_GLUT: umol_per_s {pub: in};
//var J_G_Blood: umol_per_s {pub: in};
//var v_cell: m3 {pub: in};
//var time: second {pub: in};-->
<!-- //ode(glucose_s, time) = 1e-9{m3_mole_per_litre_per_umol}/v_cell*(-J_GLUT-J_G_Blood);
enddef;-->
<!-- def comp Blood_concentration as
var glucose_s: M {pub: in};
var glucose_B: M {init: 0.012, pub: out};
var J_G_Blood: umol_per_s {pub: out};
var K: dimensionless {init: 1e-6, pub: out};-->
<!-- //var v_cell: m3 {pub: in};-->
<!-- var time: second {pub: in};-->
<!-- J_G_Blood = K*(glucose_s-glucose_B);
ode(glucose_B, time) = 0{dimensionless};
enddef;-->
<connection>
<map_components component_1="rate_constants" component_2="phenomonological_constants"/>
<map_variables variable_1="ks_12" variable_2="ks_12"/>
<!--vars k_12 and k_12;-->
<map_variables variable_1="k_21" variable_2="k_21"/>
<!--vars k_23 and k_23;-->
<map_variables variable_1="k_32" variable_2="k_32"/>
<map_variables variable_1="k_25" variable_2="k_25"/>
<map_variables variable_1="k_52" variable_2="k_52"/>
<map_variables variable_1="k_34" variable_2="k_34"/>
<map_variables variable_1="k_43" variable_2="k_43"/>
<map_variables variable_1="k_45" variable_2="k_45"/>
<map_variables variable_1="k_54" variable_2="k_54"/>
<map_variables variable_1="k_56" variable_2="k_56"/>
<map_variables variable_1="k_65" variable_2="k_65"/>
<map_variables variable_1="k_61" variable_2="k_61"/>
<map_variables variable_1="k_16" variable_2="k_16"/>
</connection>
<connection>
<map_components component_1="phenomonological_constants" component_2="Apical_concentrations"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
<map_variables variable_1="glucose_m" variable_2="glucose_m"/>
</connection>
<connection>
<map_components component_1="phenomonological_constants" component_2="Cell_concentration"/>
<map_variables variable_1="J_Na_SGLT" variable_2="J_Na_SGLT"/>
<map_variables variable_1="J_Gl_SGLT" variable_2="J_Gl_SGLT"/>
</connection>
<connection>
<map_components component_1="phenomonological_constants" component_2="parameters"/>
<map_variables variable_1="k0_23" variable_2="k0_23"/>
<!--vars C_T and C_T;-->
<map_variables variable_1="F" variable_2="F"/>
</connection>
<connection>
<map_components component_1="phenomonological_constants" component_2="Apical_voltage"/>
<map_variables variable_1="I_Na_SGLT" variable_2="I_Na_SGLT"/>
<!--vars I_Gl_SGLT and I_Gl_SGLT;-->
</connection>
<connection>
<map_components component_1="rate_constants" component_2="Apical_concentrations"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
<map_variables variable_1="glucose_m" variable_2="glucose_m"/>
</connection>
<connection>
<map_components component_1="rate_constants" component_2="Cell_concentration"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="glucose_i" variable_2="glucose_i"/>
</connection>
<connection>
<map_components component_1="rate_constants" component_2="parameters"/>
<map_variables variable_1="k0_12" variable_2="k0_12"/>
<map_variables variable_1="k0_21" variable_2="k0_21"/>
<!--vars k0_23 and k0_23;-->
<map_variables variable_1="k0_32" variable_2="k0_32"/>
<map_variables variable_1="k0_25" variable_2="k0_25"/>
<map_variables variable_1="k0_52" variable_2="k0_52"/>
<map_variables variable_1="k0_34" variable_2="k0_34"/>
<map_variables variable_1="k0_43" variable_2="k0_43"/>
<map_variables variable_1="k0_45" variable_2="k0_45"/>
<map_variables variable_1="k0_54" variable_2="k0_54"/>
<map_variables variable_1="k0_56" variable_2="k0_56"/>
<map_variables variable_1="k0_65" variable_2="k0_65"/>
<map_variables variable_1="k0_61" variable_2="k0_61"/>
<map_variables variable_1="k0_16" variable_2="k0_16"/>
<map_variables variable_1="delta" variable_2="delta"/>
<map_variables variable_1="alpha_p" variable_2="alpha_p"/>
<!--vars alpha_pp and alpha_pp;-->
<!--vars C_T and C_T;-->
<map_variables variable_1="n" variable_2="n"/>
<map_variables variable_1="z_c" variable_2="z_c"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="mu" variable_2="mu"/>
</connection>
<connection>
<map_components component_1="CFTR" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Cl" variable_2="z_Cl"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="CLC2" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Cl" variable_2="z_Cl"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="CFTR" component_2="Cell_concentration"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="HCO3_i" variable_2="HCO3_i"/>
<map_variables variable_1="HCO3_m" variable_2="HCO3_m"/>
<map_variables variable_1="J_CFTR" variable_2="J_CFTR"/>
<map_variables variable_1="J_CFTR_HCO3" variable_2="J_CFTR_HCO3"/>
</connection>
<connection>
<map_components component_1="CLC2" component_2="Cell_concentration"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<!-- vars HCO3_i and HCO3_i;
vars HCO3_s and HCO3_s;-->
<map_variables variable_1="J_CLC2" variable_2="J_CLC2"/>
<!--vars J_CFTR_HCO3 and J_CFTR_HCO3;-->
</connection>
<connection>
<map_components component_1="CFTR" component_2="Apical_concentrations"/>
<map_variables variable_1="Cl_m" variable_2="Cl_m"/>
</connection>
<connection>
<map_components component_1="CLC2" component_2="Basolateral_concentrations"/>
<map_variables variable_1="Cl_s" variable_2="Cl_s"/>
</connection>
<connection>
<map_components component_1="CFTR" component_2="Apical_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
<map_variables variable_1="I_CFTR" variable_2="I_CFTR"/>
<!--vars I_CFTR_HCO3 and I_CFTR_HCO3;-->
</connection>
<connection>
<map_components component_1="CLC2" component_2="Basol_voltage"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
<map_variables variable_1="I_CLC2" variable_2="I_CLC2"/>
<!--vars I_CFTR_HCO3 and I_CFTR_HCO3;-->
</connection>
<connection>
<map_components component_1="ENaC" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="ENaC" component_2="Cell_concentration"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="J_ENaC" variable_2="J_ENaC"/>
</connection>
<connection>
<map_components component_1="ENaC" component_2="Apical_concentrations"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
</connection>
<connection>
<map_components component_1="ENaC" component_2="Apical_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
<map_variables variable_1="I_ENaC" variable_2="I_ENaC"/>
</connection>
<connection>
<map_components component_1="IK" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="IK" component_2="Cell_concentration"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="J_IK" variable_2="J_IK"/>
</connection>
<connection>
<map_components component_1="IK" component_2="Basolateral_concentrations"/>
<map_variables variable_1="K_s" variable_2="K_s"/>
</connection>
<connection>
<map_components component_1="IK" component_2="Basol_voltage"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
<map_variables variable_1="I_IK" variable_2="I_IK"/>
</connection>
<connection>
<map_components component_1="BK" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="BK" component_2="Cell_concentration"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="Ca_i" variable_2="Ca_i"/>
<map_variables variable_1="J_BK" variable_2="J_BK"/>
</connection>
<connection>
<map_components component_1="BK" component_2="Apical_concentrations"/>
<map_variables variable_1="K_m" variable_2="K_m"/>
</connection>
<connection>
<map_components component_1="BK" component_2="Apical_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
<map_variables variable_1="I_BK" variable_2="I_BK"/>
</connection>
<connection>
<map_components component_1="NaKCC" component_2="Cell_concentration"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="J_NaKCC" variable_2="J_NaKCC"/>
</connection>
<connection>
<map_components component_1="NaKCC" component_2="parameters"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="GLUT2"/>
<map_variables variable_1="glucose_i" variable_2="glucose_i"/>
<map_variables variable_1="J_GLUT" variable_2="J_GLUT"/>
</connection>
<connection>
<map_components component_1="Basolateral_concentrations" component_2="GLUT2"/>
<map_variables variable_1="glucose_s" variable_2="glucose_s"/>
<!--vars J_GLUT and J_GLUT;-->
</connection>
<!-- def map between interstitium_concentration and parameters for
vars time and time;
vars v_cell and v_cell;
enddef;-->
<connection>
<map_components component_1="parameters" component_2="GLUT2"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<!-- def map between GLUT2 and Basol_voltage for
vars I_GLUT and I_GLUT;
enddef;-->
<!-- def map between Cell_concentration and Na_Cl for
vars Na_i and Na_i;
vars Cl_i and Cl_i;
vars J_NaCl and J_NaCl;
enddef;-->
<connection>
<map_components component_1="Cell_concentration" component_2="NHE3"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="H_int" variable_2="H_int"/>
<map_variables variable_1="H_Ext" variable_2="H_Ext"/>
<!-- vars HCO3_i and HCO3_i;
vars HCO3_m and HCO3_m;-->
<!-- vars Cl_i and Cl_i;-->
<map_variables variable_1="J_NHE3_Na" variable_2="J_NHE3_Na"/>
<map_variables variable_1="J_NHE3_H" variable_2="J_NHE3_H"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="NaK_channel"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="NaK_channel"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="J_Na_NaK" variable_2="J_Na_NaK"/>
<map_variables variable_1="J_K_NaK" variable_2="J_K_NaK"/>
</connection>
<connection>
<map_components component_1="Basolateral_concentrations" component_2="NaK_channel"/>
<map_variables variable_1="K_s" variable_2="K_s"/>
</connection>
<connection>
<map_components component_1="Basol_voltage" component_2="NaK_channel"/>
<map_variables variable_1="I_Na_NaK" variable_2="I_Na_NaK"/>
<map_variables variable_1="I_K_NaK" variable_2="I_K_NaK"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="NBC"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<!-- vars Cl_i and Cl_i;-->
<map_variables variable_1="J_NBC_Na" variable_2="J_NBC_Na"/>
<map_variables variable_1="J_NBC_HCO3" variable_2="J_NBC_HCO3"/>
<map_variables variable_1="HCO3_i" variable_2="HCO3_i"/>
<map_variables variable_1="HCO3_s" variable_2="HCO3_s"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="NBC"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<connection>
<map_components component_1="Basol_voltage" component_2="NBC"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
<map_variables variable_1="I_NBC_HCO3" variable_2="I_NBC_HCO3"/>
<map_variables variable_1="I_NBC_Na" variable_2="I_NBC_Na"/>
</connection>
<connection>
<map_components component_1="Basolateral_concentrations" component_2="NBC"/>
<map_variables variable_1="Na_s" variable_2="Na_s"/>
</connection>
<!-- def map between Cell_concentration and NBC for
vars Na_i and Na_i;
enddef;-->
<connection>
<map_components component_1="Cell_concentration" component_2="Apical_Diffusion"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="J_D_Na" variable_2="J_D_Na"/>
<map_variables variable_1="J_D_K" variable_2="J_D_K"/>
<map_variables variable_1="J_D_Cl" variable_2="J_D_Cl"/>
</connection>
<!-- def map between Apical_concentrations and Cell_concentration for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Apical_concentrations" component_2="Apical_Diffusion"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
<map_variables variable_1="K_m" variable_2="K_m"/>
<map_variables variable_1="Cl_m" variable_2="Cl_m"/>
<!--vars F and F;-->
</connection>
<connection>
<map_components component_1="parameters" component_2="Apical_Diffusion"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="z_Cl" variable_2="z_Cl"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="P_Na" variable_2="P_Na"/>
<map_variables variable_1="P_k" variable_2="P_k"/>
<map_variables variable_1="P_Cl" variable_2="P_Cl"/>
</connection>
<!-- def map between parameters and Cell_concentration for
vars A_Apical and A_Apical;
enddef;-->
<!-- def map between parameters and Apical_voltage for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Apical_Diffusion" component_2="Apical_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="Apical_voltage"/>
<map_variables variable_1="capacitance" variable_2="capacitance"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="Basolateral_Diffusion"/>
<map_variables variable_1="Na_i" variable_2="Na_i"/>
<map_variables variable_1="K_i" variable_2="K_i"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="J_D_Bl_Na" variable_2="J_D_Bl_Na"/>
<map_variables variable_1="J_D_Bl_K" variable_2="J_D_Bl_K"/>
<map_variables variable_1="J_D_Bl_Cl" variable_2="J_D_Bl_Cl"/>
</connection>
<!-- def map between Basolateral_concentrations and Cell_concentration for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Basolateral_concentrations" component_2="Basolateral_Diffusion"/>
<map_variables variable_1="Na_s" variable_2="Na_s"/>
<map_variables variable_1="K_s" variable_2="K_s"/>
<map_variables variable_1="Cl_s" variable_2="Cl_s"/>
<!--vars F and F;-->
</connection>
<connection>
<map_components component_1="parameters" component_2="Basolateral_Diffusion"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="z_Cl" variable_2="z_Cl"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Basol" variable_2="A_Basol"/>
<map_variables variable_1="P_Bl_Na" variable_2="P_Bl_Na"/>
<map_variables variable_1="P_Bl_k" variable_2="P_Bl_k"/>
<map_variables variable_1="P_Bl_Cl" variable_2="P_Bl_Cl"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="Cell_concentration"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
<map_variables variable_1="time" variable_2="time"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="NHE3"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
<!-- vars time and time;
vars R and R;
vars T and T;-->
</connection>
<connection>
<map_components component_1="Apical_voltage" component_2="NHE3"/>
<!--vars v_sc and v_sc;-->
<map_variables variable_1="I_NHE3_H" variable_2="I_NHE3_H"/>
<map_variables variable_1="I_NHE3_Na" variable_2="I_NHE3_Na"/>
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="NHE3"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
</connection>
<!-- def map between Cell_concentration and NBC for
vars Na_i and Na_i;
enddef;-->
<!-- def map between AE1 and NBC for
vars HCO3_s and HCO3_s;
vars HCO3_i and HCO3_i;
enddef; -->
<connection>
<map_components component_1="parameters" component_2="AE1"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
<!--vars A_Apical and A_Apical;-->
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="AE1"/>
<map_variables variable_1="Cl_i" variable_2="Cl_i"/>
<map_variables variable_1="HCO3_i" variable_2="HCO3_i"/>
<map_variables variable_1="HCO3_m" variable_2="HCO3_m"/>
<map_variables variable_1="J_AE1_Cl" variable_2="J_AE1_Cl"/>
<map_variables variable_1="J_AE1_HCO3" variable_2="J_AE1_HCO3"/>
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="AE1"/>
<map_variables variable_1="Cl_m" variable_2="Cl_m"/>
</connection>
<!-- def map between Apical_voltage and AE1 for
//vars v_sc and v_sc;-->
<!-- vars I_AE1_Cl and I_AE1_Cl;
vars I_AE1_HCO3 and I_AE1_HCO3;
enddef;-->
<!-- def map between Basolateral_concentrations and Basol_voltage for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Basolateral_Diffusion" component_2="Basol_voltage"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
</connection>
<connection>
<map_components component_1="parameters" component_2="Basol_voltage"/>
<map_variables variable_1="capacitance" variable_2="capacitance"/>
<map_variables variable_1="time" variable_2="time"/>
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="Na_m" variable_2="Na_m"/>
<map_variables variable_1="K_m" variable_2="K_m"/>
<map_variables variable_1="Cl_m" variable_2="Cl_m"/>
<!--vars F and F;-->
</connection>
<!-- def map between Apical_concentrations and parameters for
vars time and time;-->
<!-- //vars F and F;
enddef;-->
<connection>
<map_components component_1="parameters" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="F" variable_2="F"/>
<map_variables variable_1="z_Na" variable_2="z_Na"/>
<map_variables variable_1="z_k" variable_2="z_k"/>
<map_variables variable_1="z_Cl" variable_2="z_Cl"/>
<map_variables variable_1="R" variable_2="R"/>
<map_variables variable_1="T" variable_2="T"/>
<map_variables variable_1="A_Apical" variable_2="A_Apical"/>
<map_variables variable_1="P_P_Na" variable_2="P_P_Na"/>
<map_variables variable_1="P_P_k" variable_2="P_P_k"/>
<map_variables variable_1="P_P_Cl" variable_2="P_P_Cl"/>
</connection>
<connection>
<map_components component_1="Cell_concentration" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="J_D_P_Na" variable_2="J_D_P_Na"/>
</connection>
<!-- def map between Apical_concentrations and Paracellular_voltage for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Paracellular_Diffusion" component_2="Paracellular_voltage"/>
<!-- vars I_D_P_Na and I_D_P_Na;
vars I_D_P_K and I_D_P_K;
vars I_D_P_Cl and I_D_P_Cl;-->
<map_variables variable_1="v_ms" variable_2="v_ms"/>
</connection>
<!-- def map between parameters and Paracellular_voltage for
vars capacitance and capacitance;
vars time and time;
enddef;-->
<connection>
<map_components component_1="Basol_voltage" component_2="Paracellular_voltage"/>
<map_variables variable_1="v_sc" variable_2="v_sc"/>
</connection>
<connection>
<map_components component_1="Apical_voltage" component_2="Paracellular_voltage"/>
<map_variables variable_1="v_mc" variable_2="v_mc"/>
</connection>
<connection>
<map_components component_1="Apical_voltage" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="I_D_P_Na" variable_2="I_D_P_Na"/>
<map_variables variable_1="I_D_P_K" variable_2="I_D_P_K"/>
<map_variables variable_1="I_D_P_Cl" variable_2="I_D_P_Cl"/>
</connection>
<connection>
<map_components component_1="Basol_voltage" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="I_D_P_Na" variable_2="I_D_P_Na"/>
<map_variables variable_1="I_D_P_K" variable_2="I_D_P_K"/>
<map_variables variable_1="I_D_P_Cl" variable_2="I_D_P_Cl"/>
</connection>
<connection>
<map_components component_1="Basolateral_concentrations" component_2="Paracellular_Diffusion"/>
<map_variables variable_1="Na_s" variable_2="Na_s"/>
<map_variables variable_1="K_s" variable_2="K_s"/>
<map_variables variable_1="Cl_s" variable_2="Cl_s"/>
<!--vars F and F;-->
</connection>
<!-- def map between interstitium_concentration and Blood_concentration for
//vars J_G_Blood and J_G_Blood;-->
<!-- vars glucose_s and glucose_s;
enddef;-->
<!-- def map between parameters and Blood_concentration for
vars time and time;
enddef;-->
<connection>
<map_components component_1="Cell_concentration" component_2="A_GLUT2"/>
<map_variables variable_1="glucose_i" variable_2="glucose_i"/>
<map_variables variable_1="J_A_GLUT" variable_2="J_A_GLUT"/>
</connection>
<connection>
<map_components component_1="Apical_concentrations" component_2="A_GLUT2"/>
<map_variables variable_1="glucose_m" variable_2="glucose_m"/>
<!--vars J_GLUT and J_GLUT;-->
</connection>
<connection>
<map_components component_1="parameters" component_2="A_GLUT2"/>
<map_variables variable_1="v_cell" variable_2="v_cell"/>
</connection>
<!-- def map between Cell_concentration and Hydration for
vars H_int and H_int;
vars HCO3_i and HCO3_i;
vars J_buf_C and J_buf_C;
vars H_Ext and H_Ext;
vars HCO3_m and HCO3_m;
vars J_buf_A and J_buf_A;
enddef;-->
<!-- def map between parameters and Hydration for
vars time and time;
vars v_cell and v_cell;
enddef;-->
<!-- def map between Paracellular_Diffusion and Hydration for
vars CO2_i and CO2_i;
vars CO2_m and CO2_m;
vars J_D_P_CO2 and J_D_P_CO2;
enddef;-->
</model>
